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A recent theory for the ordered phase of hehcal or chiral magnets such as MnSi is used to cal- 
culate observable consequences of the helical Goldstone modes or helimagnons. In systems with no 
quenched disorder, the helimagnon contributions to the specific heat coefficient is shown to have 
a linear temperature dependence, while the quasi-particle inelastic scattering rate is anisotropic in 
momentum space and depends on the electronic dispersion relation. For cubic lattices the generic 
temperature dependence is given by a non-Fermi-liquid T'^/^ behavior. The contribution to the 
temperature dependence of the resistivity is shown to be r^/^ in a Boltzmann approximation. The 
helimagnon thus leads to nonanalytic corrections to Fermi-liquid behavior. Implications for experi- 
ments, and for transport theories beyond the Boltzmann level, are discussed. 
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I. INTRODUCTION 

The itinerant helical magnet MnSi has attracted a con- 
siderable amount of interest lately. This material shows, 
at ambient pressure, helical magnetic order below a criti- 
cal temperature Tc ~ 30 K J. The wavelength of the helix 
is 2tt /q, with q « 0.035 A the pitch wave number ji Appli- 
cation of hydrostatic pressure p monotonically decreases 
Tc until Tc vanishes at p = Pc ~ 14kbar.^ A tricritical 
point is observed on the phase boundary at T « 10 K, 
such that the paramagnetic-to-helimagnetic transition at 
higher pressures and lower temperatures is of first order, 
while the transition at lower pressures and higher temper- 
atures is of second or very weakly first order i^i^ In the or- 
dered phase, neutron scattering shows the helical pattern 
of the magnetization with the axis of the helix pinned 
in the [lll]-direction due to crystal field effects^ In the 
disordered phase, quasi-static remnants of helical order 
are still observed at low temperatures close to the phase 
boundary. The distribution of the helical axis orienta- 
tion is much more isotropic than in the ordered phase, 
with broad maxima in the [110]-direction.^ Such rem- 
nants of order in the disordered phase are not entirely un- 
expected close to a first order phase transition boundary.— 
In the entire disordered phase, up to a temperature of a 
few Kelvin, and up to the highest pressures investigated 
(« 2pc), pronounced non- Fermi- liquid behavior of the 
resistivity is observed, with the temperature dependence 
of the resistivity given by p{T ^ 0) ~ const. + T^/"^^ 
In the ordered phase, on the other hand, the transport 
is observed to be Fermi-liquid-like, with a leading T^- 
dependence of the resistivity. No explanation has been 
given so far for the unusual properties in the disordered 
phase. However, it is natural to speculate that the rem- 
nants of the helical order that are clearly observed in the 
paramagnetic phase have something to do with them. As 



a first step in investigating this possibility, the effects of 
the helical magnetic order on the electronic properties in 
the ordered phase of an itinerant electron system need to 
be understood. 

In a previous paper, hereafter denoted by I, we pre- 
sented a theory for the long-ranged order and fluctua- 
tions in the helically ordered phase of itinerant chiral 
magnets.— In particular, we obtained the helical Gold- 
stone modes or helimagnons in the ordered phase. Phys- 
ical quantities computed included the leading contribu- 
tion to the dynamical magnetic susceptibility at wave 
numbers near the helical pitch wave number, and the 
noninteracting single-particle Green function in the or- 
dered phase. Results for the helimagnon consistent with 
I were obtained in Ref. I^- 

In the present paper we use the results of I to deter- 
mine some of the thermodynamic and transport proper- 
ties of the helical magnetic state. In particular we will 
calculate the effects of the magnetic fluctuations in the 
ordered phase on the specific heat coefficient and on the 
electrical resistivity to first order in these fluctuations. 
We will see that the magnetic soft modes lead to non- 
analytic corrections to the standard Fermi-liquid theory 
results. Specifically, we find that the specific heat co- 
efficient 7 has a contribution linear in the temperature 
T, whereas Fermi-liquid theory gives a leading correc- 
tion to the constant Pauli value that is proportional to 
T^ InT. Similarly, the resistivity we find to have a T^/^ 
term while Fermi-liquid theory gives a T^ correction to 
the leading T^ contribution. Finally, the quasi-particle 
or inelastic relaxation rate has a temperature dependence 
proportional to 7^3/2^ which is stronger than the T^ be- 
havior found in Fermi-liquid theory. Similar effects in 
metallic antiferromagnets, in particular an anomalously 
large scattering rate, have been discussed in Ref. 

The plan of this paper is as follows. In Section II 



we give simple physical plausibility arguments that show 
how our results arise and how they are related to the 
behavior of more familiar systems. In Section III the 
results are derived using a combination of field theory, 
many-body perturbation theory, and transport theory. 
In Section IV we further discuss our results and their 
experimental implications. 



II. SIMPLE PHYSICAL ARGUMENTS, AND 
RESULTS 

In this Section we use simple physical arguments to 
obtain some of our results and discuss certain aspects of 
them. More detailed and thorough technical derivations 
will be given in later sections. 



A. Specific lieat 

In I we showed that in the chiral magnetic state there 
is one Goldstone mode, the helimagnon. Denoting the 
frequency of the helimagnon by loq and the wave vector 
by k, the dispersion relation in the long- wavelength limit 
is given byiS 



tribution. From Eq. (|2.2p we obtain 



ujQ{k) 



Zzkl + CiA;1 



(2.1) 



where Cz and c_l are elastic constants which can be ex- 
pressed in terms of the exchange splitting, the pitch wave 
number, and the Fermi wave number. The wave vector 
k = {kj_,kz) has been separated into longitudinal and 
transverse components with respect to the pitch wave 
vector q — {0,0, q), which we take to point in the z- 
direction. Notice that this dispersion relation is strongly 
anisotropic, and softer in the direction transverse to the 
pitch vector, ujq ^ k\ for k^ — 0, than in the z-direction, 
loq ^ kz for k^ = \k±\=0. 

To the extent that the helimagnon is a well-defined 
quasi-particle, one expects its contribution to the internal 
energy density to be 



u[T) 



^^wo(p)ri(u;o(p)). (2.2) 



Here n{x) = l/(exp(x/r) — 1) is the Bose distribution 
function. V denotes the system volume, and throughout 
this paper we use units such that k^ = h = 1. By using 
Eq. (12.11) and scaling out the temperature it is easy to see 
that the leading helimagnon contribution to the specific 
heat, C = du/dT, is proportional to T^jii Adding the 
leading Fermi liquid contribution, which is linear in T, 
the low-temperature specific heat in the helical phase is 



C(r^O) = 7r + 72 T2 + C'(T^ In T). 



(2.3) 



Here 7 is the usual Fermi- liquid specific heat coefficient, 
and 72 is the coefficient of the leading helimagnon con- 



1 



72 



3C(3) 



(2.4) 



with C(x) denoting the Ricmann zeta-function. The term 
of O(r^lnr) in Eq. (|2.3p is the usual nonanalytic term 
in Fermi-liquid theory that also exists in a nonmagnetic 
metal.- ^ 

It is also interesting to compare the helimagnon con- 
tribution to the specific heat contribution from the Gold- 
stone modes or spin waves in antiferromagnets and fcrro- 
magnets, which have dispersion relations wo(p) oc |p | 
and ujo{p) (X p^ , respectively. Equation (|2.2p yields 
Cafm oc T^ and Cfm oc T^/^ for these two cases. ^^ The 
helimagnetic case is thus in between the ferromagnetic 
and antiferromagnetic cases, as one would expect based 
on the nature of the respective Goldstone modes. 



B. Quasi-particle relaxation rate, and electrical 
resistivity 

The other main result of the present paper concerns the 
temperature dependence of the quasi-particle relaxation 
rate and the electrical resistivity due to the scattering of 
electrons by helimagnons in a helically ordered magnetic 
state. 



1. Quasi-particle relaxation rate 

To make plausible our result for the quasi-particle re- 
laxation rate, let us recall the case of electrons with an 
energy-momentum relation efc that are scattered by an 
effective dynamical potential V(jj, ifln), with p a mo- 
mentum and Qn = ^nTn a bosonic Matsubara frequency. 
The quasi-particle relaxation rate is given by the imag- 
inary part of the self energy, and to lowest order in the 
potential the latter is given by the diagram shown in Fig. 
[1] A standard calculation yields the following relaxation 



-^- 



FIG. 1: The lowest-order electronic self energy E in terms of 
the potential V (dashed line) and the electronic Green func- 
tion (solid line). 

rate for a quasi-particle with wave vector k on the Fermi 
surface, 



1 

Tfc 



-E 



V"ip-k,^p) 



sinh(^p/r) 



(2.5) 



where V"{p,uj) — ImV {p,iiln -^ ijJ- 
of the dynamical potential, and ^p = 
Fermi energy. 



' iO) is the spectrum 
£p ^ £f with e-p the 



In the case of a Fermi liquid, the relevant effective po- 
tential is the dynamically screened Coulomb interaction, 
which has the property V"{p,uj) oc ui/pvp!^ The inte- 
gration along the Fermi surface then gives simply a ge- 
ometrical prefactor, and the temperature dependence of 
the relaxation rate is given by the integration over the 
modulus of p,^^ 



— oc / d^ 



? 



sinh(e/T) 



(xT^ 



(2.6) 



Another well-known example is the scattering of elec- 
trons by acoustic phonons, in which case V"{p,ll!) oc 
cp5{u> — cp)r^ with c the speed of sound. This leads 
to 



1 



^c-ph 



^E 



V 
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We note that the delta-function contribution to the ef- 
fective potential reflects the phonon propagator, whereas 
the prefactor cp reflects a matrix element squared that 
describes the coupling between electrons and phonons. ^^ 
Also, we have ignored a second contribution to the 
phonon propagator that leads to the same temperature 
dependence. 

With these preliminary considerations one can make an 
educated guess about the helimagnon contribution to the 
quasi-particle scattering rate. First of all, the quantity 
^p, which measures the distance from the Fermi surface, 
needs to be generalized to reflect the existence of two 
non-spherical Fermi surfaces split by a Stoner gap. It 
follows from I Eqs. (4.13), and will be discussed in more 
detail in Sec. IIIII below, that this generalization reads 



Cp -^ ^iAp) = 2 



Cp + Cp+q ± V (Cp ^ Cp+q)^ 



4A2 



(2.8) 

Here we neglect spin-orbit coupling effects which can 
qualitatively modify the band structureJ^ A is the ex- 
change splitting or Stoner gap from I Eq. (4.10c). 
The Fermi surfaces are defined as the loci of wave vec- 
tors k with (jji.2(fc) = 0. The inelastic lifetimes for 
quasi-particles on either Fermi surface have the same 
temperature-dependence, and for definiteness we will 
concentrate on the Fermi surface given by uji{k) — 0. 
The effective potential will be proportional to the spin 
susceptibility, which in turn is known from I to be propor- 
tional to the Goldstone propagator. For power-counting 
purposes, the spectrum of the latter is adequately repre- 
sented by (see I Eqs. (4.33)) 



p,cj) oc 



1 



^o(p) 



5{uj-uja{p)) 



(2.9) 



Here 4> is the order parameter phase defined in I (see Eqs. 
(3.4a), (4.17), (4.28), and (4.32b) in I). In addition, we 
expect the effective potential to contain a multiplicative 
function g{k,p) that describes the coupling between elec- 
trons and helimagnons, in analogy to the phonon case. 



Altogether we expect for the inverse lifetime of a quasi- 
particle with wave vector k on the 1-Fermi surface due 
to scattering by helimagnons 



1 



1 
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g{k,k + p) 



V ^ ujo{p)Bm\i{ua{p)/T) 



S{uji{k+p)-ujo{p)). 
(2.10) 



The determination of the function g requires some infor- 
mation about the coupling mechanism, as in the electron- 
phonon case. Since is a phase mode, one expects the 
physically relevant correlation function to describe the 
fluctuations of the gradient of (p, rather than of the phase 
itself. This suggests g{k,p) oc (fe — p)^. Indeed, we will 
see in Sec. IIIII that for free electrons, e^ = k'^/2mc, one 
has g{k, k + p) oc p^ for small deviations p from the 
Fermi surface. The p]_ term is absent in this case, as 
it is in the helimagnon dispersion relation. At the same 
time, cji(fc+p) = kj_-p±/mc. Power counting then shows 
that l/rfc oc T^/^ for generic points on the Fermi surface. 
However, for a generic e^ that is consistent with a cu- 
bic symmetry (as is relevant for, e.g., MnSi) one finds 
g{k,k + p) (X {k± ■ p_l)^ -I- 0{k^). At the same time, 
uji{k + p) is modified such that the leading argument of 
the delta- function is no longer proportional to k± ■ p± . 
Since, in a scaling sense, p]_ ^ Pz ^ T^ this reduces the 
power of temperature by one, and for generic points on 
the Fermi surface one obtains 



1 



'^c-hni 



2^3/2 



(2.11) 



for the electron-helimagnon scattering contribution to 
the quasi-particle scattering rate. Notice that this is 
stronger than the standard Fermi-liquid contribution to 
the scattering rate, Eq. (|2.6p . but not strong enough to 
destroy the Fermi liquid. 

Equation (|2.10[) also yields a qualitative result for 
antiferromagnetS(^ whose Goldstone modes have an 
isotropic dispersion relation with coq (p) oc |p | . Power 
counting with g{k,p) oc (A; — p)^ yields 1/re-afm oc T^. 
The temperature dependence is weaker than in the result 
for the helimagnet, as one would expect from the nature 
of the respective Goldstone modes. 



2. Electrical conductivity 

The quasi-particle lifetime is related to, but not the 
same as, the relevant time scale for the electrical con- 
ductivity a, or the resisitivity p — 1/a. Technically, the 
conductivity is given by a four-fermion correlation func- 
tion, and vertex corrections enter in addition to the self- 
energy contributions that determine the quasi-particle 
lifetime. Physically, backscattering events contribute 
more strongly to the resistivity than forward scattering 



events. This leads to a transport scattering rate l/rtr oc p 
that is given by Eq. (^31) . or the first hne in Eq. (PTU]) . 
with an extra factor of {p — k)^ in the integrand. In ad- 
dition, the transport rate gets averaged over the Fermi 
surface. In the Coulomb case this changes just the ge- 
ometric prefactor, but not the temperature dependence. 
The Fermi-liquid contribution to the resistivity thus is^ 



III. EFFECTS OF HELIMAGNONS ON 
OBSERVABLES 

In this section we use the field theory of fluctuations 
in the helically ordered phase developed in I to calculate 
the specific heat and the quasi-particle relaxation time. 



(xT\ 



(2.12) 



A. Specific heat 



In the case of scattering by a propagating mode, as in the 
phonon and helimagnon cases, the momentum is slaved 
to the frequency by the delta-function, and hence the 
temperature dependence does change. In the phonon 
case, where the phonon wave number scales as T, one 
obtains the familiar Bloch-Griineisen result^^ 



Po-ph oc T^ 



(2.13) 



In the helimagnon case, the smallest (in a scaling sense) 
additional factor in the integrand is proportional to (p — 
fc)^ ~ T, and the averaging over the Fermi surface is 
not important since the behavior at generic points on 
the Fermi surface is the leading one. We thus expect 



Pc-hr 



(xT 



5/2 



(2.14) 



Adding the Fermi-liquid contribution, we thus obtain for 
the low-temperature behavior of the resistivity in the or- 
dered phase of a clean (no impurity scattering) helimag- 
net 



p{T ^ 0) = P2 T2 + P5/2 T5/2 + 0{T^ 



(2.15) 



with p2 and p^/2 temperature-independent coefficients. 

For antiferromagnets the additional factor of (p — fc)^ 
in the integral leads to an additional factor of T^ in 
the resistivity. The contributions from the Goldstone 
modes to the resistivity in this case is thus expected to 
be Pafm oc T^^ Again, the behavior is weaker than in 
the helimagnetic case. 

The above simple arguments capture the structure of 
the full theoretical development in Sec. lIIIi although they 
ignore many subtleties that occur in the actual calcula- 
tion. The result shows that helimagnon scattering, while 
leading to a nonanalytic temperature dependence to the 
resistivity, is weaker than the usual electron-electron con- 
tribution. Also, the T^/"^ contribution is one power of T 
weaker than the observed temperature dependence in the 
paramagnetic phase of MnSii^ Nevertheless, it is intrigu- 
ing that it is a nonanalytic term and involves a half- 
integer power of the temperature. It is also intriguing 
that the temperature dependence of the quasi-particle 
rate in the helimagnetic phase is the same as the observed 
transport rate in the paramagnetic phase. We will come 
back to the possible relevance of these observations for 
the paramagnetic phase in Sec. IIVI 



In this subsection we calculate the contribution to the 
free energy from Gaussian helimagnon fluctuations and 
use the result to obtain the leading fluctuation contribu- 
tion to the specific heat. We will reproduce the results 
of the simple argument given in Sec. |TT1 

If we ignore all degrees of freedom other than the Gold- 
stone mode g, we have from Eqs. I (4.33) the following 
Gaussian action, 

^^'^ = ^ E E 5(P' *") r(P, i^) 9{-P, -m, (3.1a) 

p iSl 

with a vertex function, 

T{p,in) - -{tnf + c^2(p) + \n\j{p). (3.ib) 

Here, and in the remainder of this section, we suppress 
the integer index on Matsubara frequencies. uJo{p) is the 
oscillator frequency given by Eq. (12. ip . and the damping 
coefficient 7(p) is given by Eq. I (4.33c). In this Gaussian 
approximation, the Goldstone mode contribution to the 
grand canonical potential S is 



S(2) = -- In 
V 



D[g] cxp(^-A^'->[g]) 



This gives 



T 



-^'^ = ^EEi-r(p,.f^). 



2V 



(3.2) 



(3.3) 



p jO 



Neglecting the ifl = term, which does not contribute 
to the specific heat, the frequency sum is conveniently 
performed by using the identity (see Appendix E]) 



Y^ 



r> inr(p,ir!) = 



in 



-T 

■K 

xlm 



In sinh In IljI 

2T ' ' 



dij 



— Inr(p,c^ + z0)). (3.4) 



The specific heat at constant volume is obtained from ^ 
via the relation^i 



Cv 



-Td^E/dT^ 



_d_ 
df 



E-T 



dE_ 
df 



(3.5) 



In the limit of negligible damping, 7(p) — > 0, the expres- 
sion for Cv can be written, by combining Eqs. (j3.4p and 
^3, 



CviT)^-^l,J2''oiP)nB{u:oip))- 



(3.6) 



This is identical with the result for Cy obtained from 
Eq. (|2.2p . and leads to Eq. (|2.4p for the specific heat 
coefficient. With a nonzero damping coefhcient, the fre- 
quency integral cannot be done in closed form, but the 
asymptotic low-temperature behavior is still proportional 

toT2.25 

The model calculation given in I, which kept modes 
with wave numbers fc = and k — q, results in a heli- 
magnon dispersion relation (see Eq. I (4.33b) and Ref. 



^o(fc) = A -±-^kl/{2k^Y + i fci/(2gfcF)2. 



(3.7a) 



Keeping modes with higher wave numbers changes the 
dispersion relation to 



ujQ{k) = A 



g 

3fcF 



kl/{2k^Y + -ky{2qk^Y. (3.7b) 



This result, which is analogous to the one obtained for 
cholesteric liquid crystals, ^^ leads to the following values 
for the elastic constants in (|2.ip . 






(3.8a) 
(3.8b) 



Here fcp is the Fermi wave number, and the energy scale 
A is the exchange splitting or Stoner gap as defined in I. 
This results in the following asymptotic low-temperature 
behavior of the helimagnon contribution to the specific 
heat 



C\r{T)^q^Ac{TlT,f 



(3.9) 



and Tq — 



with a coefficient Ac = V6C(3)/27r « 0.47 

Xq'^/Qkp. Equation p.9p is valid for temperatures T <^ 

Tq. We will provide a discussion of this result, as well as 



a semi-quantitative analysis, in Sec. ITVl 



B. Quasi- particle relaxation time 

In this subsection we calculate the temperature depen- 
dence of the quasi-particle inelastic relaxation time in a 
helimagnet, i.e., the lifetime of a free-electron state on 
the Fermi surface due to scattering by helimagnon fluc- 
tuations. We will calculate this quantity to first order in 
the helimagnon susceptibility. We will use the result to 
determine the low-temperature behavior of the electrical 
resistivity. 



1. Effective action 

To proceed, we need an electronic action that takes into 
account the helical magnetic order and helical magnetic 



fluctuations. The starting point is an electronic action of 
the form 



S'[V',V'] = So[ip,il;] + - 



dxdyn\{x) Aij{x 



y)ni{y). 

(3.10a) 
Here V' and '0 a-re fermionic (i.e., Grassmann- valued) 
fields that depend on space, time, and a spin index 
a; nl{x) = ?/'Q(a;)tT^^'0/3(a;), where a* (i = 1,2,3) de- 
notes the Pauli matrices, are the components of the 
electronic spin density field ns{x). x = (a;,r) com- 
prises the position x in real space and the imaginary 
time r, and the related integration measure is given by 
J dx — Jy dx J^ dr. So contains all parts of the action 
other than the spin-triplet interaction. The interaction 
amplitude A consists of a point- like Hubbard interaction 
with an amplitude Ft and a chiral part with coupling 
constant c whose origin has been explained in I, 

Aij{x - y) = S{x - y) [5ij Ft -f eijfc cFt dk\ ■ (3.10b) 

Note that A is static; it depends only on spatial position, 
and not on imaginary time. 

The general idea is now to replace one of the spin den- 
sity fields in the last term in Eq. (|3.10ap by a classical 
(i.e., c-number valued) field that represents the effective 
field seen by the electrons due to the magnetic order. If 
one is just interested in incorporating static helical order, 
one can implement a mean-field approximation by replac- 
ing either one of the spin density fields in Eq. (|3.10ap by 



its average value, n^ 
effective action 



2n(n) 



This yields an 



So[ip,il^] = So[ip,ip] 



dx Ho{x) ■ ns{x), (3.11a) 



where 



Ho{x) = rt{n,{x)), 



(3.11b) 



and we have dropped a constant contribution to the ac- 
tion. The remaining question regards the equation of 
state that determines (rig (a;)). This has been given in I 
in a saddle-point approximation, which yields 



Ho{x) = A (cos(q • cc), sin(q • x), 0) , 



(3.11c) 



Here q is the pitch vector of the helix, and the ampli- 
tude A is determined by a generalized Stoner equation of 
state. So is the reference ensemble action given by Eq. I 
(4.7a), which describes noninteracting electrons (or elec- 
trons interacting via a spin-singlet interaction only) in an 
effective magnetic field given by the helically modulated 
magnetization. 

More generally, one might ask whether one can re- 
place Ho in Eq. (j3.11ap by a fluctuating classical field 
H{x) — Ft M{x), where M{x) represents the spin den- 
sity averaged over the quantum mechanical degrees of 
freedom. For the sake of presentational simplicity we will 
explain this procedure for an ordinary Hubbard interac- 
tion (c = in Eq. (|3.10bp and q = in Eq. ((XTTc)) ). The 



development for the chiral interaction is exactly analo- 
gous. 

Writing M{x) — Ho/Ft + 5M{x), and ignoring a con- 
stant contribution to the action, one has 



5'[V', V, SM] = So [^, V'] + Tt / dx 6M{x) 



(5ns (x), 

(3.12a) 
with Sns{x) — ns(x) — {ns{x)). This needs to be supple- 
mented by an action governing SM. If the expectation 
value of M is to represent the exact magnetization, this 
action must be 

A[6M] = —Jdxdy 5A'h{x) [xs'),, {x,y)SM,{y), 

(3.12b) 
where Xs is the spin susceptibility of the system. A purely 
electronic effective action is now obtained by integrating 
out the fluctuations 5M. The Gaussian integral yields 



Scs[ip,-ip] = So[ip,'ip]- 



dx dy Snl (x) Xs' {x,y) Snl (y) . 



(3.13) 
With the full chiral interaction, Eq. (|3.10bp . one obtains 
formally the same result. The only reservation is that in 
the chiral case the correlation function {SM SM) = Xs 
is, strictly speaking, not the physical spin susceptibility. 
However, its leading hydrodynamic contribution is the 
same as that of the latter, as was shown in I. The result 
(|3.13[) thus consists of the reference ensemble described 
by So, whose Green function has been determined in I, 
and an effective interaction given by the spin susceptibil- 
ity in the helically ordered phase. The latter has been 
calculated in I in a Gaussian approximation. Notice that 
the effective potential depends on imaginary time or fre- 
quency. 

It is obvious that the above heuristic considerations 
represent, in a diagrammatic language, some kind of in- 
finite resummation, which means that the resulting ef- 
fective action is valid only in conjunction with certain 
constraints. To clarify this point it is useful to rewrite 
the starting point, Eq. (|3.10ap . in the form 

S[^,il;]=So[^,^] + ^JdxdySnl{x)A,^ix~y)Sni{y). 

(3.14) 
Now consider the bare spin-triplet interaction A, Fig. 
[21 and its renormalization by the ladder resummation 
shown in Fig. [3] The result of this resummation is 



>--< 



G 



FIG. 2: Diagrammatic representation of the bare spin-triplet 
interaction A and the reference ensemble Green function G. 



an effective interaction of the form given in Eq. (j3.13p . 
with a Gaussian approximation for Xs- The constraint 
that must be used in conjunction with the effective action 




FIG. 3: Effective spin-triplet interaction resulting from a lad- 
der resummation. 



()3.13p is now obvious: The effective interaction must not 
be used in ways that constitute renormalizations of Xs, as 
doing so would result in double counting. In particular, it 
is safe to use Eq. (j3.13p in any perturbative calculation 
to linear order in Xs- We will now proceed and use it 
to calculate the quasi-particle inelastic lifetime to that 
order. 



2. Electronic self energy and Green function 

We consider the self energy S of the single-particle 
Green function to first order in the perturbing potential 
Fj Xs = Xs- There are two self-energy diagrams to this 
order, namely, the direct or Hartree and the exchange or 
Fock contributions shown in Figs. Ufa) andHfb), respec- 
tively. Both self-energy contributions are nondiagonal in 




(h) 



FIG. 4: The direct (a) and exchange (b) contributions to the 
single-particle self energy. 

both spin and momentum space. It is readily seen that 
the direct contribution is purely real, and hence does not 
contribute to the scattering rate. The exchange contri- 
bution is given by, 



(3.15) 



V 






Here ioj denotes a fermionic Matsubara frequency, and 
Gkp(i^) is the single-particle Green function of the ref- 
erence ensemble which is given explicitly by Eqs. I (4.13). 
All components of the self energy can be obtained from 
Eq. (|3.15p by means of straightforward, albeit lengthy, 
calculations. 

We now consider the Dyson equation for the renormal- 
ized Green function Q, which reads 

Qkp{^^) = Glliiu) ~ Efcp(z^). (3.16) 

From Eqs. I (4.13) for G and Eq. ((XTC)) it follows that E 



has a structure very similar to that of G ^ , 

Sfcp(«a;) = dkp [(7+-'>^++{k,iuj) + a-+T. — {k,iuj)] 
+ <>k+q,p 0"+ S+_ (fc, iuj) 
+ 6k-q.p 0-- S-+(fc, iuj). (3-17) 

Here the notation is the same as in Eqs. I (4.13). In 
particular, a± = (cti ± jct2)/2, a^ — — cr+cr-, and a — i- — 
a^a+, with ai^2 Pauli matrices. From this expression it 
follows in turn that Q has the same structure as G, 

Gkpii^) = 6kp[<7+- A+{k,iuj) + <7-+A-{k,iuj)] 
+ Sk+q,p cr+ B+{k, iuj) 
+ 6k-q,p(T- B-{k,iuj). (3.18a) 



Here 



A±{k,iuj) = --^ 
B±{k, iuj) 



g '^{k±q,iuj) 



g^^{k ± q, iuj) gj.^{k, iuj) — (A±(fc, iu)) 

(3.18b) 
— A±(fc,iw) 



9^^{k ± q, iuj) g±^{k, iuj) - (A±(A;, iu)f 

(3.18c) 



with 



5±\A;, iw) = G"i(fc, icj) - S±±(fc, iw), (3.18d) 



A±(fc, iw) = A — S±zp(fc, iu 



(3.18e) 



These expressions constitute an exact inverse of Eq. 
p.l6p . as can easily be checked by a direct multiplica- 
tion. 



3. Quasi-particle relaxation time 

The quasi-particle relaxation time is determined by the 
imaginary parts of the poles of the Green function G, Eqs. 
(I3.18P . For a vanishing self energy there are poles at 



UJ 



1,2 



(fe) 



a+a±q±v(^fc-^fc±'j) +4A2 



with ^fe as defined after Eq. (|2.5p . The poles with differ- 
ent signs of the square root reflect the Stoner splitting of 
the Fermi surface into two sheets. On a given sheet, lu^ 



and w are related by uj^ {k + q) 



'{k) 



1,2). AU 



poles can thus be expressed in terms of 



wi,2(fc) = 2 ( ^'^ 



Cfe+q ± V (^fc ^ ^fe+q) 



AX' 



(3.19) 
Both of these resonance frequencies are real, reflecting 
the fact that the quasi-particles are infinitely long lived to 
zeroth order in the effective potential Xs- To first order in 
Xs the resonance frequencies acquire an imaginary part, 
corresponding to a finite relaxation time T(fc), in addition 
to a shift of the real part. For definiteness, we consider 
the resonance at uji{k); the relaxation time on the other 
sheet has the same temperature dependence. We find 



r(fc) 



Im <^ i;++(fc, z) + E__(fe + q, z) 



2A 



(Cfe - ik+q) 



{ik - ^k+qf + 4A2 



YH [S++(fc,z)-S__(fc + q,^)] 



(a-a+q) +4A2 



— [I]+_(fc,z) + E_+(fc + q,z)] 



(3.20) 



r 



Here the frequency z is given by the 'on-shell condition' 

Z — LL!l{k) -\- lO. 

Of particular interest is the relaxation rate for quasi- 
particles on the Fermi surface, which is defined by 



wi(fe) 



2 ( Cfc + ^k+q 



(Cfc-a+q) +4A2 



= 0. 



1/2 



(3.21a) 



For A = q = this reduces to the usual definition of the 
Fermi surface, S,k = 0. In the general case it implies in 
particular 

aa+q^A^. (3.21b) 



With this condition, we use Eq. (|3.15p in Eq. (|3.20p . and 
the results of I for Xs, keeping only the leading hydro- 
dynamic contributions to the latter. Performing the fre- 
quency summations we find for the relaxation rate of a 
quasi-particle on the Fermi surface 



-2A2 



a:^ 



9ik,p) 



Mk) (a + ^k+qf V ^ sinh (a;i(p)/T) 

xx'l^{p-k,LUiip)). (3.22a) 

This expression is valid to determine the leading low- 
temperature behavior of 1/tf only. We have anticipated 
the fact that the dominant contribution to the latter 



comes from momenta p « A;, and accordingly have re- 
placed p by A; in all contributions to the integrand that 
have a finite and nonzero limit as p -^ k. Here Xaa 
is the spectrum of the susceptibility of the phase fluc- 
tuation variable (p that was introduced in I. Its leading 
hydrodynamic part is also proportional to the Goldstone 
mode. From I (4.33a) or, alternatively, from I (4.29a) 
and (4.40) r^ we find for the phase susceptibility 



{p.m 



1 



1 



2iVF 3fc|, uj^{p) - (ifiy 



(3.22b) 



where ujq{p) is the hclimagnon resonance frequency, see 
Eqs. (Pl^ and I (4.33b). The spectrum is thus 



1 



(P,^) - :t^^T7T 



2A^F 3fc| 2^0 (p) 



[S {uj - ujoip)) 
-5{iu + ujo{p))]- (3.22c) 



The function g is given by 



g{k,p) = uJi{p) iCk + ^k+q) 



h ^k {S,k+q — (,p+q) 

-a+q(a-^p) • (3.22d) 



Notice that g{k,k) — 0, that x"{Pi^) is soft at p = 
and u! — LUi{k) — 0, and that Eq. (|3.22ap has indeed the 
form of the educated guess, Eq. (|2.10p . 

For later reference we also introduce the 'off-shell' rate 
l/T{k, e) obtained by putting z = e + iO in the self en- 
ergies in Eq. p.20p . With k on either Fermi surface one 
finds 



-27rA2 



nik,e) {Ck + ^k+^ 



■I^Hu) + ne + u)]l,J2 



Here Wi(p) is given by Eq. p.l9p . and 

g,{k,p; e) = (w,(p) - ^p) (C^+q - e) 



9t{k,p;e)x'l^{p~k,u)6{u + e-~uj^{p)), {i ^ 1,2) , 

(3.23a) 



(^.(P) - Cp+q) i^^ + 2^ a + q 

e - t,k+q 



(3.23b) 



is a generalization of Eq. (j3.22dp . with the property 
gi(fc,p;e = 0) =g{k,p). 

To evaluate the final integral we consider the limit 
ep 3> A ^ qwF, as we did in I. The result will obvi- 
ously depend on the direction of fc. It also depends on 
the functional form electronic energy momentum relation 

Efc- 

a. Isotropic energy-momentum relation Let us first 
consider a nearly-free electron approximation, with 



Efc = k'^/2me, 



(3.24) 



with Too the effective mass of the electrons. Then 
g{k,k-\-p) (xpl -\-0{pI), and cji(fc+p) = fc^ ■p^/mc + 
0{hp,,kl). 

For fc^ = the leading result as T -^ is 



.^VeKi 



rp2 



T{k) ml kl{^u+ik+qY 



ex A 



* eF\2 f T"- ^ 



with Tq the temperature scale defined after Eq. (|3.9p . For 



fcj^ 7^ the asymptotic result is 



r(fc) 



Ar 



oc A 



Al/2qfc3 fcp 



2^5/2 



kl_ {ik + £.k+q 



m 



5/2 



(3.25b) 



with At « 289. In the second lines in Eqs. (|3.25ap and 
(j3.25bp we have assumed a generic value of k, and have 
neglected all numerical prefactors as well as terms small 
in qvp/X for clarity. The crossover between these two 
types of behavior occurs at a temperature proportional to 
Tx = X {q/kF)^{k±/kzY. It is easy to check that adding 
an isotropic quartic term, proportional to (k'^Y, to e^ 
does not change the powers of the temperature in these 
results. 

6. Cubic energy-momentum relation In an actual 
metal, the underlying lattice structure causes the elec- 
tronic energy-momentum relation to be anisotropic. In 
the case of a cubic lattice, as in MnSi, any terms consis- 
tent with cubic symmetry are allowed. For instance, to 
quartic order in k the following function is allowed, 



ffc = k^ /2mc 



'' (pj2 



'•c'^F 



'^y'^z 



^z '^x ) 7 



(3.26) 



with V a dimensionless measure of deviations from a 
nearly-free electron model. Generically, once expects 
V = 0(1). Other quartic terms that are consistent with 
a cubic symmetry, e.g., the cubic anisotropy k^ + k^ + k'^, 
can be obtained by adding an isotropic (fe^)^ term to Eq. 
(|3.26[) . For v ^ the asymptotic behavior of g is changed 
compared to the nearly-free electron case. 



g{k,k+p) = - 



f 2v 



\ ZiTThQfvT^ 



— {ki_ ■ p^) 



(Cfc + Cfc+q) 

and for q in (0,0, l)-dircction one has 

uJi{k+p) == (fc^ ■P±)/Mo 



(3.27) 



2i^ 



ZillliQlXiT^ 



2v 



ZiTi1j(\I\iT^ 



{kl + kl + k,^q + q^/2)k.^p^ 



{kl + kl + k,q + q^l2)kypy 



where 



l/^c = l/We + 1^ o„. 7 2 



o(pi,p.), 

q{2kz + q) f fc - ^k+q 



2mofcp ^fc + (k+q 

If we define 

and 

Biy^ {2k, + qf q^k^til 



(3.28a) 



(3.28b) 



V^lc 



rrickp ^ 



{kl, + kl + k,q + q^l2), (3.29a) 



Cfc = 

with 
B = 



A^ 



1,5 



Kp 



(a+Cfc+q)^ 



(3.29b) 



48V3 



dxdz 



21/433/4 y^ ^-^ ^"- Vz2 + a;4 sinh V z^ + x^ 
= 54.99... (3.29c) 

the evaluation of Eqs. (|3.22p now yields 



r(fc) 



— Ck 



{klAl + klA^fn \-X 



3/2 



oc t^ A 



if}' 



3/2 



(3.29d) 



In a cubic system, the quasi-particle relaxation rate 
thus generically shows a T^/'^ behavior, except on four 
special lines in fc-space, k^ = 0, ky = 0, and k^ = iky, 
where the prefactor of the T'^/^ vanishes and the behavior 
is T^. These results are valid for q in (0,0, l)-direction. 
The corresponding expressions for q in (1, 1, l)-direction, 
which is relevant for MnSi,^* can be obtained by a corre- 
sponding rotation of fe. That is, {k^, ky,kz) in the above 
expressions should be replaced by 



kx -^ [kx - VSky + V2kz)/V6, 
ky -^ {kx + viky + v2fc2)/v6, 
k^ -^ i-2kx + V2k,)/V6. 



(3.30) 



It needs to be stressed that the T^/^ and T"^/^ terms 
due to the helimagnons, as well as the Fermi-liquid T^ 
term, all contribute to the quasi-particle relaxation rate, 
and which of these contributions dominate in a given 
temperature regime is a quantitative question. The same 
is true for the different contributions to the specific heat. 
In Sec. IIV Al we will give estimates for parameter values 
that are appropriate for MnSi. 



C. Resistivity 

The electrical resistivity can be obtained as the inverse 
of the conductivity, which in turn is given by the Kubo 
formula 



cTijiiQ) = — — [Trij{ifl) — TTij{ifl — 0)] . (3.31a) 



Here 



7T,,{in) - -e^Tj2 ^E^»('^)"^-(P) 



V 

ni,n2 k.p 

X (V'ni,o-(fe)V'ni+«,cr(fc)V'n2,(T'(p)V'n2-n,o-'(p)) • 

(3.31b) 

is the current-current susceptibility tensor or polarization 
function, with v{k) = dck/dk. The average, denoted by 
(...), is to be performed with the effective action given 
in Eq. (|3.13p . The simplest approximation to tt is a fac- 
torization of the four-point correlation function in Eq. 
()3.29bp into two Green functions. The difference between 
this simple approximation and the full polarization func- 
tion is customarily expressed in terms of a vector vertex 
function T with components ^F, 



Ti^j{in) = -e^T^ — ^^,F^^(iw,itJ-if7) 
X E E *^^('^') S^k'i^^) Sk^pi^^ - ^^), (3.32) 



fe' o- 



with Q from Eq. p.l6p . It is well known that care must 
be taken to use consistent approximations for the self en- 
ergy S that defines Q and the vertex function. ^^-"^^ The 
simplest combination that fulfills the consistency require- 
ment, which is equivalent to the Boltzmann equation for 
the conductivity, is a self-consistent Born approximation 
for the self energy, 



^fpi^^) 



^E^EEC''(-+^^) 

«0 k'p' a' [3' 



'XVaa',l3'l3{k' ~ k,p' -pjil), (3.33a) 
and a ladder approximation for the vertex function. 
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Ti^{tuj,iuj-in) = SkpS^ptv{k) + Tj2^ Yl Ti,',^',',{iuj + in',tLj-tn + in')g^''^!{iuj + in') 



V 

iO,' k'. k", p', p" 

Qp'i' (*^ -i^ + i^') I [V^'c M' {k k',p p '; i^') + Vf3p> „.„(p ' - P , fe' - fc; i^')] . (3.33b) 



r 



Here V is the potential in the effective action, 

Vaf3,jsik, p;in) = Y a'^p x'J (fe, p ; ifl) a^^g, (3.33c) 



with Xs the spin susceptibility times T^ as used in the 
calculation of the quasi-particle relaxation time in Sec. 
IIIIB| see Eq. p.lSp . Equations (|3.3ip are represented 
diagrammatically in Fig. \5\ With the results for the 



0- 





< 



FIG. 5: Consistent approximation for the self energy E and 
the vertex function F. 

spin susceptibility given in I (4.41) - (4.43) we find that 
only four matrix elements of the potential V are nonzero. 
They can be expressed in terms of the phase susceptibil- 
ity, Eq. (l3:22bl) . 

Vi2, 21 {k,p;i^) = 6kpX'^X4>4>{k-q,i^), (3.34a) 
^21,21 (fe,p;^f^) = SkpX^xM'« + q,i^), (3.34b) 

^12,12(fc,P;i^) = -'5fc-2q,pA2x00(fc-q,i^), 

(3.34c) 

V21,2l(fe,P;i^) = ~Sk+2q.pX'^X<p<pik + q,'i^), 

(3.34d) 

The development of the transport theory now proceeds 
analogously to the electron-phonon scattering case, with 



playing the role of the dynamical potential. The chief 
complications are, (1) the spin structure, which leads to 
two independent matrix elements, T^^ and F^^ of the 
vertex function, and (2) the fact that the vector vertex 
function Tkp depends on the pitch vector q in addition to 
k. {p gets eliminated by the Kronecker-(5 constraints in 
Eqs. (|3.18a[) and (|3.34[) .) In addition, there are two Fermi 
surfaces, a feature the helimagnetic problem shares with 
the ferromagnetic one. 

As a consequence of the dependence of the vertex func- 
tion on q, the polarization function, and hence the con- 
ductivity, is no longer isotropic. The tensor is still di- 
agonal, but the transverse and longitudinal components, 
with respect to q, are different, 

TTy (if7) = Si-j [{5ii -f 5i2) TTTiiVl) + (5i37rL(iri)] (3.36) 

As a consequence of the two Fermi surfaces, the conduc- 
tivity is the sum of two contributions. 



- Ji) 



.(2) 



CTL.T = Re lim crL,T(«f^ -^uj + iQ) = aY-'+a]^'. (3.37) 

After calculations that are quite involved, but follow the 
same reasoning as in the electron-phonon case, one ob- 
tains the following expression for the transverse conduc- 
tivity. 



^i) 



2mi 



Jde[-^)^T.'(Mk))kl 



XTj(fc,e)Aj(/c,e). 



(3.38a) 



X'^X4>4>iP^^^) = yiP,i^) 



(3.35) 



Here r^ is the quasi-particle relaxation rate given by Eqs. 
(I3.23p . LJi are the resonance frequencies given by Eq. 
(1^ . and /(e,T) = l/exp(e/T-|-l) is the Fermi distribu- 
tion function. The scalar vertex functions A^ are related 
to the part of the vector vertex function T^p that is pro- 
portional to k. As in the electron-phonon case, they obey 
integral equations which reduce to algebraic equations for 
the purpose of determining the leading temperature de- 
pendence only of the conductivity. We obtain 



A,(fc,e) = l-A,(fc,6) J^^'^^''\, /'duK«) + /(e + M)]1^5.(fc,P;e)/(fe,P)^"(fc-p,«)<5(w-t^,(p)), (3.38b) 



where 



/(fc,p) = l-^^(fc-p) 



(3.38c) 



Substituting the solution of Eq. (|3.38bp into Eq. p.38ap 
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shows that the conductivity is given in terms of trans- 
port relaxation time Ttr that is different from the quasi- 
particle relaxation time r. Namely, 



i^) eW f-df\ 1 



with 



-5]<5(c.,(fc))fcir«(fe,6), 
(3.39a) 



-27r 



2r,V^^(fc,e) i^+g+^kyJ ^ 



du 



— [n{u) + /(6 + u)] 



17 X! 9t{k,p) V^'ik - p, u) 6{u - UJ^{p)). 

(3.39b) 



V 



Here 



V^'{k^p,u) 



jk-p) 
2kl 



V"{k 



p,u) 



(3.39c) 



From Eqs. (f3?34l) and (I3.22cp it follows that {k ~ p)l^ - 
M ~ T in a scaling sense. The leading temperature de- 
pendence of the transport relaxation rate 1/rtr will there- 
fore carry one higher power of the temperature than the 
quasi-particle relaxation rate 1/t. 

The leading contribution to the longitudinal conduc- 
tivity is given by Eq. (|3.38ap . with k'^ replaced by k^. 
Since k is pinned to the Fermi surface by the (5-function 
in Eq. ()3.38a|) . this difference does not change the temper- 
ature dependence. The leading temperature dependence 
of (7l is thus also again given by rtr . 

Combining these results and observations with the re- 
sults of Sec lIIIlTSl wc find the following asymptotic tem- 
perature dependence for both the transverse and longi- 
tudinal transport relaxation rates due to helimagnons in 
the ordered phase of a helimagnet with cubic symmetry: 



l/4cxl/r-(x -^ r-^ 



T 



T„ 



5/2 



T 



T„ 



7/2' 



kp) \\ 

(3.40a) 
Here v is the prefactor of the cubic anisotropy in Eq. 
(|3.26p . Tq is defined after Eq. p.9p . and we have again 
omitted all numerical prefactors. For both the longitu- 
dinal and transverse resistivities, this leads to Eq. (|2.15p 
with 



P5/2 oc i^"* A 



kv) \\) 



1 



t; 



,5/2 



(3.40b) 



IV. DISCUSSION AND CONCLUSION 

In summary, we have calculated the effects of the Gold- 
stone mode in the ordered phase of helical magnets, the 
helimagnon, on the low-temperature behavior of the spe- 
cific heat, the quasi-particle relaxation rate, and the re- 
sistivity. The detailed microscopic calculations given in 



Sec. IHII have corroborated the simple physical plausibil- 
ity arguments given in Sec.|TTl The helimagnon contribu- 
tion to the specific heat was found to have a T^ temper- 
ature dependence, whereas the corresponding contribu- 
tion from ferromagnetic and antiferromagnetic Goldstone 
modes goes as T^/"^ and T"^, respectivelyj^ The quasi- 
particle relaxation rate was found to have a T^/"^ tem- 
perature dependence, which is nonanalytic and stronger 
than the T^ behavior in a Fermi liquid. The resistiv- 
ity depends on the transport relaxation rate, whose tem- 
perature dependence is one power weaker than that of 
the quasi-particle rate. The asymptotic low-temperature 
contribution to the resistivity due to electron- helimagnon 
scattering thus goes as T^/^. 

We divide our discussion of these results into a semi- 
quantitative discussion that makes predictions of experi- 
mental relevance, and more general theoretical remarks. 



A. Predictions for Experiments 

The leading low-temperature results given in Sec. IIIII 
hold for temperatures T <^ Tq, with Tq the temperature 
scale given after Eq. (|3.9p . This can be seen as follows. 
The helimagnon dispersion relation as given in Eq. (j2.ip , 
or (|3.7bp . is valid for wave numbers k < q, and crosses 
over to a different behavior around k = q. The energy or 
temperature scale related to this crossover is thus given 
by 



Tq =uja{k^ ^ q,k±=0) = ^Xq^ /k^ 



(4.1a) 



Here we have defined a dimensionless number C2 by c^ — 

Vithin ou: 
we have C2 = 1/36, or 



C2 X^q^ /k^. Within our weak-coupling model calculation 



Tq^Xq^/&kl, 



(4.1b) 



which is also the energy of a ferromagnon at fc = g within 
Stoner theory. Estimates for the parameters entering 
Tq that are appropriate for MnSi have been given in I: 
g/fcp « 0.024, and ep ~ 23, 000 K. The value of A is less 
certain, for reasons explained in I. An upper limit is pro- 
vided by the Fermi energy, and we use A/ep = 1/2 as a 
practical upper bound. A likely lower limit is given by 
A — 520 Kj^ Given the large ordered moment in MnSi)^ 
a value of A that is a sizeable fraction of ep is actually 
more plausible than the latter. For these two extreme 
choices one obtains Tq w 1.1 K and Tq k, 50 mK, respec- 
tively. In providing such estimates one should keep in 
mind that even semi-quantitative estimates are difficult 
for a system like MnSi that is characterized by a com- 
plex Fermi surface and strong interactions.'^^ It should 
also be noted that Tq gives only the general scale for the 
crossover; a numerical evaluation of, for instance, the in- 
tegral that determines the specific heat, Eq. p.6p . shows 
clear deviations from the T^ behavior already at a tem- 
perature of about Q.2Tq. The relevant temperature scale 
for the interesting helimagnon effects is thus quite low. 
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To see what the behavior crosses over to at T « Tg , we 
reahze that the term k\ in Eq. (j3.7bp is the leading (in a 
scahng sense) contribution to a term k^. For wave num- 
bers k larger than q, the dispersion relation is given by 
uJo{k) ex k^, and fcz and k± both scale effectively as T^/^. 
For the specific heat at T > Tg this implies a helimagnon 
contribution proportional to T^^^. This is the same be- 
havior as in a ferromagnet, since the dispersion relation is 
effectively ferromagnet-like and the specific heat depends 
only on the dispersion relation. 

A scaling expression for the helimagnon contribution 
to the specific heat that incorporates both sides of this 
crossover is 



both helimagnon contributions to the relaxation rate is 



C^^{T)^q^fciT/T,). 



(4.2) 



Here fc is a universal scaling function with fc{x — » 0) = 
Ac x^ and fci^ ^ 1) = Be x^l"^ . Here Aq is given after 
Eq. (j3.9p . and Be is another universal number. 

It is also interesting to compare the helimagnon contri- 
bution to the Fermi-liquid one. In the case of the specific 
heat, this means comparing the first term in Eq. (|2.3p 
with the second one. In a nearly-free electron model, the 
coefficient 7 of the Fermi-liquid term in Eq. (|2.3p is given 
by 7 = fc^/eep.- From Eq. ([SJ]) . we have 72 = q^Ac/T"^. 
The helimagnon contribution is equal to the Fcrmi-liquid 
one at a temperature 



T* = 7/72 



1 



J- a 



36yl(7 ep 9 



(4.3) 



With the value of k-p/q appropriate for MnSi given above, 
and A a sizeable fraction of ep, one finds T* Ri Tg; for 
smaller values of A, T* is correspondingly smaller than 
Tq (which itself has a smaller value) . Notice that T <T* 
does not preclude an observation of the T^ helimagnon 
contribution to the specific heat; it can be extracted by 
plotting C/T versus T. 

For the quasi-particle relaxation rate, which can be 
measured by means of either weak-localization or tun- 
nelling experiments, a similar discussion applies. The 
results given in Sec. IIIIB 3l applv for temperatures small 
compared to Tg, and at higher temperatures a crossover 
occurs to a regime that is characterized by an effectively 
quadratic helimagnon dispersion relation and k^ ^ k± ^^ 
T^/^. Repeating the analysis in Sec. IIII B 31 shows that 
in this regime, 1/t(A;) cx T. (We note in passing that 
this is not the correct result for a ferromagnet, where the 
Goldstone dispersion relation is effectively the same, but 
the structure of the equation analogous to Eq. (|3.22ap 
is different.) This is true irrespective of whether one 
uses the nearly-free electron energy-momentum relation 
of Sec. lIHBSal of the cubic one of Sec. IIII B 3 bl This 
linear temperature dependence is remarkable, although 
this is not the asymptotic low-temperature contribution, 
since it deviates so strongly from the quadratic Fermi- 
liquid result. 

A scaling expression that incorporates this limit for 



r{k) 



(4.4) 

with /^^/^^(x ^ 0) (X a;^/2, f^^^'^\x -^ 0) cx x^/^, and 

/^^/'^(x > 1) cx /^^/'^(x > 1) cx x. 

Again, it is instructive to compare with the corre- 
sponding Fermi-liquid result, which is 



T 
rpL 8 V ep 



(4.5) 



We start with the T^/^ ^gj.j^^ ^^ (|3.25bp . which will be 
present even in systems with a sizeable lattice anisotropy. 
It is larger than the Fermi-liquid contribution for temper- 
atures larger than a crossover temperature 



T 



3^6 



2-5/2 - 2AE 



T 



(4.6) 



Due to the high power of the uncertain parameter A/ep 
that enters Eq. (|4.6p it is hard to give even a semi- 
quantitative estimate. For A a sizeable fraction of ep this 
result suggests that T^_r^j2 is substantially larger than 
Tg, which means that the T^/^ contribution, in the re- 
gion of its validity, is small compared to the Fermi-liquid 
term. However, for system with a smaller value of A/ep 
the T^/^ contribution may be important. We now turn 
to the T^/'^ law that is the asymptotic helimagnon con- 
tribution to the relaxation rate for a cubic system. By 
comparing Eqs. (|3.29p with Eq. (|4.5p one finds that the 
T"^/^ law becomes comparable to the Fermi-liquid T^ be- 
havior at a temperature 



T 



3/2-2 



6 



d>Bv^ 



kv) V A J 



T„ 



(4.7) 



0.025 Tg. Smaller values of A 
'3/2-2/^9- ^3/2-2 ^ -^9 i^ neces- 



with B from Eq. (|3.29cp . The T^^^ term is numerically 
large compared to the Fermi liquid T^ contribution for 
temperatures T < T^/2-2- "^^^ numerical value of Tg*,2_2 
depends both on a high power of the parameter v in addi- 
tion to the high power of ep/A. With ly = 1 and ep/A = 2 
as above one finds T'3*/2_2 
result in a larger ratio T 

sary for the T^''^ term to be numerically large compared 
to the Fermi-liquid contribution in the entire range of its 
validity, and for the change from the asymptotic T^/^ 
behavior to the pre-asymptotic T term to occur at a 
temperature low enough for the helimagnon contribution 
to be still large compared to the Fermi-liquid contribu- 
tion. According to the above estimate this is the case for 

At this point we need to remember that the helimagnon 
dispersion relation we have used to calculate the relax- 
ation rate, Eq. (|2.ip . is valid only for rotationally invari- 
ant systems, while the T^/^ behavior is the consequence 
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of an anisotropic electron dispersion relation, Eq. (|3.26p . 
However, since the modification of the helimagnon dis- 
persion is small, namely, on the order of the spin-orbit 
coupling squared, there still is a large temperature range 
where the T^/'^ behavior is realized. To see this, we recall 
that the spin-orbit coupling gso, which is proportional to 
q, changes Eq. dHH) to (see I Eq. (2.23)) 



Wo(fe) 



^fc^ + c^fcl +bczq'^k'\ 



Ikl 



(4.8) 



where 6 is a number on the order of unity. If we scale kz 
and k±_ with appropriate powers of T, Cz, and c_l, such 
that the main contribution to the integral in Eq. (|3.22ap 
comes from fc^ « /c^ ~ li we see that the spin-orbit term 
in 0^0 is unimportant for temperatures 



T > T,o - bXiq/k^f cc (q/kpfT,. 



(4.9) 



On the other hand, the T'^/^ term dominates over the 
2i5/2 ^Qj-j^ fQj. temperatures 



T < t: 



-'3/2-2^2-5/2 



%Bv^ 



kp. 



'3/2-5/2 y ^■i/2-2'2-b/2 A 

(4.10) 
As expected, Tso is small compared to T^/2-5/2 ^y ^ ^^'-^~ 
tor of [q/kpf ex g^^. 

Finally, the helimagnon contribution to the resistivity 
is suppressed compared to the contribution to the quasi- 
particle scattering rate by a factor of T/X. For realistic 
temperatures one thus expects the Fermi-liquid T^ con- 
tribution to dominate. This is consistent with the state- 
ment in Ref. y that in the ordered phase of MnSi, the 
resistivity shows a T^ behavior. 



B. General remarks 

All of our results hinge on the Goldstone mode not 
being overdamped. In an ultraclean system, a generic 
band structure actually leads to a weak overdamp- 
ing of the mode in certain directions in momentum 
space, which can change the above power laws at very 
low temperaturesi^ However, even small amounts of 
quenched disorder qualitatively weaken the damping, see 
Sec. IV. E. 2 and Ref. 34 in I. The results also hinge on ro- 
tational invariance. If the direction of the helix is pinned 
by the underlying lattice, then the y'k^+k^ dispersion 
relation crosses over to a linear one at asymptotically 
small wave numbers, and the temperature dependences 
of the observables will cross over to the same power laws 
as in the acoustic phonon case. The crossover temper- 
ature for this effect to become relevant depends on the 
strength of the pinning. 

For the specific heat we have found that the asymp- 
totic low-temperature dependence is in between the one 
for ferromagnets and antiferromagnets, respectively, as 
one might expect from the fact that the helimagnon dis- 
persion relation, Eq. (j2.ip . is antiferromagnet-like in the 



longitudinal direction, and ferromagnet-like in the trans- 
verse one. We point out, however, that the Goldstone 
susceptibility at zero frequency diverges more strongly 
with vanishing wave vector k than in the ferromagnetic 
case. Namely, from Eq. (j3.23bp (or from I Eq. (4.33a)) 
we see that in the helimagnetic case the static suscepti- 
bility diverges as l/(fcz + k\), whereas in a ferromagnet 
the corresponding behavior is 1/fc^. One expects this to 
have dramatic effects for the hydrodynamics, as is the 
case in certain liquid crystals. '^^ This problem will be 
studied separately. 

The transport theory presented in this paper is still at 
a rather primitive level. We have calculated the resis- 
tivity in the simplest possible consistent approximation, 
which corresponds to a solution of a Boltzmann equation 
where it is assumed that the bosonic modes remain in 
thermal equilibrium. Whether this zero-loop approxima- 
tion yields indeed the leading low-temperature behavior 
of the resistivity is an open question. In ferromagnets, 
mode-mode coupling effects (or loops in a field-theoretic 
language) have been found to induce a square-root fre- 
quency dependence in the zero-temperature resistivity^^ 
Possible similar effects in the helimagnetic case, and the 
nature of the mode-mode coupling effects at nonzero tem- 
perature, remain to be explored. The results of the 
current calculation, namely, effects of the helimagnons 
on the conductivity that are weaker than the usual T^ 
Fermi-liquid behavior, are consistent with reports of a T^ 
behavior of the resistivity in the ordered phase of MnSi;^ 
However, transport in the ordered phase has not been in- 
vestigated systematically, and the observation of possible 
mode-mode coupling effects may require measurements 
at very low temperatures. 

The quasi-particle relaxation rate 1/t has the some- 
what surprising property that its temperature depen- 
dence depends on the details of the electronic energy- 
momentum relation e^. As was shown in Sec. IIII B 31 
the generic behavior for an e^ consistent with a cubic 
symmetry is 1/r oc T"^'^, whereas a nearly- free electron 
model leads to t c>c T^/^ for generic k. This dependence 
of asymptotic low-temperature properties on microscopic 
details is unusual and raises questions about universal- 
ity. We also note that the prefactor of the leading T^/^ 
behavior depends on the fourth power of the parameter 
V that describes deviations from a nearly-free electron 
model, see Eq. (I3.29cp . In systems where the deviations 
from nearly-free electron behavior are small, the asymp- 
totic T'^/^ law could then be confined to extremely low 
temperatures. 

We stress that the current theory deals with the or- 
dered phase of a helimagnet, whereas experimentally 
much more interesting behavior, namely, a T"^' ^ behavior 
of the resistivity in violation of Fermi-liquid theory, has 
been observed in the disordered phases- While the cur- 
rent theory has nothing to say about this, the observed 
remnants of helical order in at least part of the region 
where the non-Fermi-liquid behavior is observed^ makes 
an understanding of the ordered phase a likely necessary 
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prerequisite for a theoretical treatment of the disordered 
phase. 

Finally, we mention that the results for the quasi- 
particle relaxation rate and the resistivity presented 
above hold for clean systems. The presence of quenched 
disorder can drastically change the transport behavior; 
this will be investigated in a separate publicationi^ 



can be written 



F=— / dLJi^{u;)f{Lj + iO), 



where 



(Alb) 
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APPENDIX A: 



MATSUBARA FREQUENCY 

SUMS 



In this appendix we derive the identity p.4p . Consider 
an even function f(z) = f{—z) of one complex argument 
z. Let / have singularities (poles or cuts) on the real 
axis, but be analytic for Im z > 0. Let ri„ = 2iTTn 
(n integer) be the set of bosonic Matsubara frequencies. 
Using standard techniques the functional 



F = r^/(zf^„) 



(Ala) 



viz) — — coth — 

^ ^ 2 2T z 



(Ale) 



The advantage of i'{z) over the more commonly used Bose 
distribution function n{z) is that ^{z) is odd in z and has 
no pole at z = 0. On the real axis, n has the property 



^H-tA LsinhM-lnlcl 



(A2) 



Using this in Eq. (jAlbp . and integrating by parts, yields 



rp /'OC 

F=^ duj 

■^ Jo 



Insinh In |cij| 

2T ' ' 



— Im/(c^ + zO). 

(A3) 

Only the imaginary part of /(w + iO) contributes since 
the real part is an even function of lu. With / = InF this 
yields Eq. 
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